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THERMAL AND DIELECTRIC PROPERTIES OF A 
HOMOGENEOUS MOON BASED ON MICROWAVE AND 
INFRARED TEMPERATURE OBSERVATIONS 

SUMMARY 


A homogeneous lunar model was used to establish theoretically the 
relationship of the thermal properties of the moon to its dielectric properties. 
This was accomplished using lunar observational data obtained in the micro- 
wave and infrared wavelengths. Measured values of loss tangent, dielectric 
constant, and thermal diffusivity for postulated lunar materials are compared 
to values calculated from microwave and infrared observations. Lunar micro- 
wave brightness temperatures are compared to subsurface temperatures based 
on selection of a thermal parameter and a depth into the lunar interior. 

The measured values of loss tangent and dielectric constant for the 
postulated materials agree with the values based on lunar observations; 
however, it is not possible to establish which sample offers the best comparison. 
The microwave brightness temperature curves do not offer good agreement with 
the calculated subsurface temperatures since the phase lag and amplitudes could 
not be made to agree simultaneously. As a result, it may be possible to use 
the homogeneous model to relate the electrical properties of the lunar material 
to its thermal properties. 


INTRODUCTION 


This report relates the dielectric and thermal properties of a homo- 
geneous moon based on calculations using microwave and infrared observational 
data of the moon. The calculations are compared to measured dielectric and 
thermal properties of materials thought to be similar to that of the moon to 
determine if the homogeneous model is sufficient for making calculations con- 
cerning the moon. 

The experimental data used in these calculations are published data of 
experimentalists working in areas of microwave and infrared astronomy, and 
thermal and dielectric properties of materials. 



LUNAR MODEL AND FOURIER HEAT EQUATION 


General Parameters 


The moon is assumed to be a smooth homogeneous solid which conducts 
heat in the manner of an isotropic body. Mathematically, the moon is con- 
sidered to be a semi-infinite solid. The surface temperature is a harmonic 
function of time and the period of rotation is full moon to full moon (2. 55 x 10 6 
seconds) . 

The model is of density p (g/cm 8 ), thermal conductivity K (W/cm °K) , 
and heat capacity C( j/g °K) and is of dielectric material with negligible mag- 
netic loss. The magnetic permeability is considered to be that of free space. 


Fourier Heat Equation 


The thermal behavior of the lunar model is represented by the Fourier 
conduction equation for one dimension in an isotropic body; hence, 


8 2 T 


1 9T 


— o - — = 0 

8x z a 8t 


(x s 0, t > 0) 


(1) 


The boundary conditions to be applied are 


(a) T is finite as x -* <*>; 

(b) T (x,t)|= T ( 0,t) = f(t) 


x=0 


where f(t) = A 0 + £ A n cos + B n sin • ) 

n=l ' \ ' 
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Using the solution of Ingersoll [ 1J , assuming the solution to be periodic, and 
applying boundary condition ( a) , the solutions to the differential equation are 
as follows: 


T t = Aje *^[cos (cot - xn/co/2o! ] 


( 2 ) 


T 2 = A 2 e x '^ aj / 2o! [sin (cot - x\W 2a;J 


(3) 


By adding the two solutions, using the principal of superposition to account 
for boundary condition (b) , the solution becomes 


T = A 0 + Ee- xW “ P cos(^ - x«/WSTp) 
n=l n \ ' 


B sin 
n 




2n7rt 

P 


W n7r/ Q! P 


) 


(4) 


The temperature is exponentially damped by e~ n7r a ^as x is increased, 
indicating lower temperatures beneath the surface. Also, the time at which 
a temperature wave at depth x + 0 reaches a maximum and the time at which 
the surface temperature reaches a maximum will be different. This lag in 
time is given by 


t 



(5) 


FOURIER SERIES APPROXIMATION TO LUNAR 
INFRARED SURFACE TEMPERATURES 


Measurements of lunar surface infrared temperatures may be repre- 
sented by 
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where A 0 , A^, and are computed for the number of terms necessary to 

represent the measured values of temperature to the accuracy desired. Calcu- 
lations of these coefficients are accomplished for discrete data using the 
following relations [ 2] : 

N 

A 0 = 1/2N £ f(t ) , (7) 

r=-N+l 


N 

A = l/N 2 f(t ) cosK t , (K=A0, N) , 
K r=-N+i r r 


N 


A n = 1/2N 2 cos N t^ , 


r=-N+l 
N 


B = l/N 2 f(t ) sinK t . 
K r=-N+l r r 


( 8 ) 


(9) 

( 10 ) 


This representation of the surface data allows calculations of temperatures 
below the surface at any depth x. An example of the approximation to the 
surface data and calculations of subsurface temperatures is shown in 
Figure 1. 


GEMERAL MICROWAVE THEORY 


Observation and analysis of the behavior of the electromagnetic radiation 
emitted from the moon enables one to determine the physical characteristics of 
the lunar material. Figure 2 [ 3] illustrates radiation from a point beneath the 
surface of the moon. As .the radiation moves from point p toward the surface, 
it is subject to both absorption and re -emission by the lunar material. Further, 
the radiation is subjected to the reflective and refractive properties of the 
lunar material. Details of these phenomena are omitted from this text and the 
reader is referred to Weaver [ 3] . 
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The intensity of radiation emitted from beneath the surface is converted 
to temperature for radio telescopes and the brightness temperature of the 
moon at any point [ 3] is represented by 


T b (v,0 ,t)= V 


R (0) + R (0 ) 
P §_ 


'T 0 + TTm cos I fit - y+ -77^ cos t 2 fit - | 2 1 + • • • Ui) 


aj(0) 


a 2 ( 0 ) 


T 0 is the constant component of the surface temperature; and T 2 are ampli- 
tudes of the first and second harmonic temperature. Additionally, a^e) and 
a2(0) are the amplitude attenuation factors with and | 2 the phase lag of the 
harmonics. 


The problem is simplified by restricting the temperature variation to 
the center of the lunar disc in which case 0 - 0. Thus, the effective brightness 
temperature becomes [ 4] : 



(1 - R x > 


A 0 + 


a t ( 0) 


cos ( fit- 


+ 


Aa 

a 2 (0) 


cos (2 fit 


U) + 



where the reflection coefficient is [ 4] 


( 12 ) 


R 


nT 7 - l 
, ^T+ l 


(13) 


and e is the relative dielectric constant of the lunar material. 
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THERMAL AND DIELECTRIC RELATIONSHI PS 


Using the solution to the Fourier heat conduction equation and the modi- 
fied equation for the intensity of radiation emitted from the moon, it is possible 
to formulate a relationship between the thermal diffusivity and the dielectric 
conductivity of a material and compute its dielectric constant and lo ss tange nt. 

The thermal wave, as shown previously, is damped by a factor e X n?r 
At some depth, x, the damping term will become e -1 . The depth x is termed 
the depth of penetration, denoted 1 , of the thermal wave and is given by the 
following expression: 


l t = (aP/ir) 2 . (14) 

Also, there is a relationship between the depth of penetration of the thermal 
wave and the depth of penetration of the electromagnetic wave; i. e. , 


1 


e 



(15) 


where 1 is the depth of penetration of the electromagnetic wave and 6 is the 
solution to the equation [ 5] 


(1 + 26 + 2 8 2 ) 2 = -Sr 1 

ie i 


Bl 1A°1 

B 0 T 0 ( 0) 


(16) 


with Te 0 and Te A being the first and second terms in a Fourier Series approxi- 
mation to the microwave temperature throughout a lunation; Bj and B e are 
averaging coefficients assumed to be unity; T 0 ( 0) and Tj( 0) are the first and 
second terms in a Fourier series approximation to the infrared surface tempera- 
tures throughout a lunation. The reflection coefficient for perpendicular incidence 
of radiation on the moon is given by 


= i T eo 
1 ' B 0 T 0 (0) * 


(17) 


6 



If the reflection coefficient for perpendicular incidence is known the dielectric 
constant may be calculated by the relation 


-(R 2 + 1) 

L<R*- i) J 


( 18 ) 


The depth of penetration of the electromagnetic wave ( attenuation distance) [ 6] 
is also given by 


1 r> 

= “ r 

e 47 r € 2 


i) "2 

■i +|[ 1 + (tan A) 2 ] 2 ) 


( 19 ) 


The solution of this equation for tan A gives the loss tangent as 


tanA= {[(X 2 /87r 2 e 1 2 ) + l] 2 - l} 2 . 

G 

In terms of the thermal diffusivity equation ( 20 ) becomes 


( 20 ) 


tan A = {[ (A. 2 / 8 tt e a 6 2 P) + 1 ] 2 - 1 } 2 . 


( 21 ) 


Thus, a relationship between the loss tangent and the thermal diffusivity is 
established and may be simplified to read 


tan A = (c^ar 2 + 2c 1 q; *) 


,-l\2 


( 22 ) 


where e* = \ 2 /87re6 2 P . 


The dielectric conductivity is given by the expression 


a = 2 ir tan A v e € 0 , ( 23 ) 

where v is the frequency of the radiation and e 0 is the dielectric constant of 
free space. 
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Substitution for tan A gives the dielectric conductivity as a function of thermal 
diffusivity, thus 

o- = K l0 T 2 + K 2 a _1 , (24) 


vW\ 


WithKl = I?$F2 and K 2 = 


'KV € 


ZLK 


62 p 


Hence, relationships are established which relate the thermal properties of a 
simulated lunar material to the dielectric properties. These relationships are 
given parametrically in Figures 3 and 4, 


The calculations for Figures 3 and 4 are based on infrared data Set I 
and Set n, respectively. Under the conditions of this model, it is possible to 
compare the behavior of the dielectric conductivity with that of thermal 
diffusivity. For example, consider the laboratory-measured value of basalt 
powder at 3.28 cm to be 1. 681 x 10 -5 cm 2 / sec. The dielectric conductivity 
then ranges from 1. 2 x 1 Or 4 to 2. 0 x 10 -4 M£2 /m when considering the mi- 
crowave measurements in the wavelength range 3. 15 to 3. 20 cm. 


It is also interesting to note that for fixed thermal diffusivity there is a 
pronounced shifting of the parameter lines toward larger values of dielectric 
conductivity as the wavelength of microwave observation decreases. Also, 
comparison of Figure 3 with Figure 4 denotes a slight shift toward larger values 
of dielectric conductivity for fixed thermal diffusivity. This shows that 
differences in surface temperatures do cause different results although 
apparently not so drastic as the changing of wavelength of microwave 
observations. 

To refine this portion of the analysis requires that selected laboratory 
samples be chosen and their thermal and dielectric properties be observed at 
predetermined wavelengths. Simultaneously, microwave observations should 
be obtained at the same wavelength. Additionally, the refinement of the surface 
temperatures throughout a lunation should be continued, especially at or near 
midnight and pre-dawn. 


EXPERIMENTAL DATA 


The experimental data used in comparisons and calculations related in 
this study are from published reports in areas of infrared and microwave 
astronomy and measurements of postulated lunar materials. 
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The infrared measurements used were published by Low [7] , Sinton [8], 
Saari [91, and Murray and Wildey [10J. These data are shown in Figures 5 and 6. 
For this study, the infrared data were combined into two sets of data. Infrared 
dataSet I consisted of the data published by Sinton, Low, and Saari. Infrared 
data Set II consisted of data published by Sinton, Low, and Murray and Wildey. 

A curve representing the surface temperatures was drawn through the data 
points and temperature values read at even intervals throughout the period. 

These values were then used in the Fourier series approximation to the lunar 
surface temperatures. Figure 7 shows the surface curve as represented by the 
Fourier series approximation to the combined data of Saari, Sinton, and Low. 

Also shown are the surface data with a difference curve between the data and the 
approximation. 

The microwave measurements used are shown in Table I [41. These 
measurements were chosen primarily because the wavelengths at which the 
observations were made are comparable to the wavelength at which the measure- 
ments of the simulated lunar materials were made. 

The A. D. Little Laboratory measurements of dielectric constant [11], 
loss tangent, and thermal diffusivity were used in order to make comparisons 
of calculated values with measured values of these quantities. The calculated 
values of dielectric constant and loss tangent were obtained using the diffusivity 
values for basalt (as measured by A. D. Little) and the microwave measure- 
ments listed in Table I. Table II shows the measured data which were compared 
to the calculated values . 


TEST OF THE HOMOGENEOUS MODEL 


The adequacy of the homogeneous model is tested according to the method 
presented by Weaver [3]. The phase lag of the first harmonic, was plotted as a 
function of the observed amplitude attenuation factor aj, where a } is given by 

1 + 2 6 + 26 2 = a 2 ! . (25) 

The calculated values should fall on the theoretical curve of Figure 8 so that the 
homogeneous model accurately represents the observations. The model is con- 
sidered inaccurate if the values of a t are consistently above the theoretical curve. 
Weaver does not consider this to be a discriminating test because of the manner 
in which the ratios A/ (Aj/aj) are established. 
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The values of phase lag are taken from the observational data. The 
factor 5 is calculated based on the observational data using equation (16) , and 
a-f is then calculated using equation (25) . These values are then plotted on the 
theoretical curve denoting the relationship of these factors. Figure 8 shows 
how the observational data compare to the theoretical values. This comparison 
denotes that the homogeneous model does fit the observational data in a repre- 
sentative manner. 


RESULTS 


Calculations of dielectric conductivity as a function of thermal diffusivity 
are shown in Figures 3 and 4. Figure 3 is based on the microwave data of the 
authors indicated on the figure, and the infrared data are shown in Figure 5. 
Figure 4 is based on the same microwave data and the infrared data of Figure 6. 

As a result of the calculations of dielectric conductivity, the equation 

_A 

or = K t a~ 2 + K^a -1 was reduced to <7 = \/ K 2 0 ! 2 , allowing a simplified relation- 

ship betwe en d ielectric conductivity and thermal diffusivity. Table III gives the 
values of n / K 2 for the microwave data considered. 


Table II compares the calculated values of dielectric constant and loss 
tangent of the homogeneous model to laboratory measurements of postulated 
lunar materials. The powders offer better comparison than the solid materials, 
but it cannot be said which of the two powders compares best for values of 
dielectric constant. Figures 9 and 10 offer a graphical comparison of loss tan- 
gent values with basalt making the better comparison. The data of Mayer , 
McCullough, and Slonaker offer the best comparison to the loss tangent values of 
basalt and pumice powders , but it is difficult to say which observational data 
best compares the calculated and measured values of dielectric constant. 

Figures 11, 12, and 13 compare the calculated and measured values of dielectric 
constant. The test of the homogeneous model indicates that the calculations of 
the amplitude attenuation factor from observational data are in general agreement 
with the theoretical curve, although some values do not appear acceptable. 

The calculations of temperature and phase lag at depths below the lunar 
surface are compared to microwave observations in Figures 14 and 15.* Each 
set of microwave observational data shown in Table III was used in turn to 
generate a brightness temperature curve throughout a lunation. Since the con- 
stant component, phase lag, and the first harmonic differ between observers, 


'"'In Figure 14, the calculated subsurface curve is based on infrared data Set I and 
is the same in all parts of Figure 14. Likewise, in Figure 15, the calculated 
subsurface curve is based on infrared Set II and is the same in all of Figure 15. 
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the generated brightness temperature curves differ accordingly. These 
brightness temperature curves were compared to the subsurface temperature 
curves calculated, based upon the infrared data of Set I and Set II, respectively. 

It was anticipated that these comparisons would allow a subsurface 
temperature curve to correspond to a brightness temperature. The brightness 
temperature curves are indicated by the observer's name and the wavelength at 
which the observations were made. The calculated subsurface temperatures . 
are indicated by type of material and the depth of penetration of the thermal 
wave. 


These comparisons are not acceptable, in general, since the microwave 
temperature values do not properly coincide with the infrared subsurface 
temperatures at a depth. The depth was selected based on solution of equation 
(14) . It should be noted that the depth of penetration of the thermal wave 
differs for basalt and pumice caused by different values of thermal diffusivity; 
however, the subsurface temperatu re curve is the same for both. This is 

because the attenuation factor, e is set equal to e -1 . 

Another unacceptable result of this comparison is that the time of 
maximum temperatures for the microwave and calculated subsurface curve do 
not coincide, or, simply, the values of phase lags do not agree. 

These results are difficult to explain since the attenuation may extend 
farther than one optical depth. The fact that the microwave temperature is an 
integrated value over depth, and the calculated surface temperature corresponds 
to a particular depth, may cause this discrepancy , or it may be that the homo- 
geneous model is inadequate. 


CONCLUSION 


The homogeneous model may be used to relate the electrical properties 
of the moon to its thermal properties. Since the surface temperatures were 
used in the initial boundary conditions, little can be said about its adequacy 
as a thermal model. The calculations of temperature and phase lag at depths 
below the lunar surface do not agree with the observed microwave temperature 
values; therefore, either the model is inadequate or the microwave data have 
not been treated properly so as to allow comparison with calculations made for 
a single point on the lunar central region. 
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The calculated values of dielectric constant and loss tangent agree in 
general with the laboratory measurements of basalt and pumice powders; how- 
ever, it cannot be definitely stated which of the samples offers the better com- 
parison. Thus, it may be impossible to determine the composition of the moon 
from dielectric properties using a homogeneous model. 

The observational microwave data yielded different values of dielectric 
constant and loss tangent for different wavelengths of observations. This may 
be due to differences in observation techniques rather than any phenomena 
associating wavelength to dielectric constant. 

It must be concluded that when treated in the manner of this analysis, 
the homogeneous model cannot be expected to generate both thermal and 
dielectric properties of the lunar material nor can it distinguish its composition 
or structure. As a result it must be considered inadequate and other models 
must be derived for these purposes. 

George C . Marshall Space Flight Center 

National Aeronautics and Space Administration 
Huntsville, Alabama, June 4, 1968 
908-20-03-0000-28-00-000 
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TABLE I. EXPERIMENTAL MICROWAVE DATA 


a 

X 

T b 

J-o 

c 

Ti 

d 

AUTHOR 

0. 13 

216 

120 

16 

Fedoseyev, 1963 

0.40 

204 

56 

23 

Kislyakov, Plechkov, 1963 

0. 80 

197 

32 

40 

Salomonovich, 1958 

0. 80 

211 

40 

30 

Salomonovich, Losovskiy, 1962 

1.63 

208 

37 

30 

Kamenskaya, Semenov, Troitskiy, 
Plechkov, 1962 

2.00 

190 

20 

40 

Salomonovich, Koschenko, 1961 

3. 15 

195 

12 

44 

Mayer, McCullough, Slonaker, 1961 

3.20 

223 

17 

45 

Losovskiy, Salomonovich, Koschenko, 
1961 

3.20 

210 

13.5 

55 

Krotikov, Porfiryev, Troitskiy, 
1961 

3.20 

213 

14 

26 

Bondar, Zelinskaya, Porfiryev, 
Strezhneva, 1962 

3.20 

216 

16 

15 

Bondar, 1962 


a. X Wavelength at which observations made (cm) 

b. T 0 Constant microwave temperature component (° K) 

c. Tj Amplitude of first harmonic (°K) 

d. £ Phase lag of amplitude of first harmonic ( deg, ) 
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TABLE IL CALCULATED AND MEASURED LUNAR MATERIAL PROPERTIES 
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indicates the wavelength at which microwave observations were made (calculations based 
on infrared data Set I) . 




TABLE II. CALCULATED AND MEASURED LUNAR MATERIAL PROPERTIES 

( Concluded) 
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TABLE IE. VALUES OF \Tk 2 FOR RELATIONSHIP BETWEEN 
DIELECTRIC CONDUCTIVITY AND THERMAL DIFFUSIVITY, 

a = \Tk 2 «~ 2 


X 

n/k 2 

AUTHOR/DATE 

0. 13 

0. 1322 x 10~ 4 

Fedoseyev, 1963 

0.40 

0.3926 X 10- 5 

Kislyakov, Plechkov, 1963 

0. 80 

0.2198 x 10" 5 

Salomonovich, 1958 

0.80 

0.2106 x 10“ 5 

Salomonovich, Losovskiy, 1962 

1.63 

0. 2058 x lO" 5 

Kamenskaya, Semenov, 
Troitskiy, Plechkov, 1962 

2. 00 

0.1466 x 10" 5 

Salomonovich, Koschenko, 1961 

3. 15 

0. 7626 x lO' 6 

Mayer, McCullough, 
Slonaker, 1961 

3.20 

0. 5164 x 10“ 6 

Losovskiy, Salomonovich, 
Koschenko, 1961 

3.20 

0. 6268 x 10~ e 

Krotikov, Porfiryev, Troitskiy, 
1961 

3.20 

0. 6044 x 10~ 6 

Bondar, Zelinskaya, 
Porfiryev, Strezhneva, 1962 

3.20 

0. 6420 x 10" 6 

Bondar, 1962 
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ALPHA OJ X 10- 3 



FIGUKE 1. FOURIER SERIES SURFACE CURVE AND 
SUBSURFACE TEMPERATURES 



FIGURE 2. MICROWAVE EMISSION FROM THE MOON 
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DIELECTRIC CONDUCTIVITY (mho/m) 



I0* 7 lO- 6 10 " 5 I0* 4 I0 -3 I0 -2 10"' 10° 

THERMAL DIFFUSIVITY (cm 2 /sec) 


NOTE: The number preceding observer's name denotes wavelength 
in centimeters. 


FIGURE 3. DIELECTRIC CONDUCTIVITY VERSUS THERMAL 
DIFFUSIVITY BASED ON INFRARED DATA SET I 
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DIELECTRIC CONDUCTIVITY { mho/m) 



NOTE: The number preceding observer’s name denotes wavelength 
in centimeters. 


FIGURE 4. DIELECTRIC CONDUCTIVITY VERSUS THERMAL 
DIFFUSIVITY BASED ON INFRARED DATA SET II 
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TIME (FRACTION OF A PERIOD) 

FIGURE 5. LUNAR SURFACE TEMPERATURE DATA OF 
LOW, SINTON, AND SAARI 



0 0.1 0.2 0.3 0.4 0.5 06 07 0.8 09 1.0 

TIME (FRACTION OF A PERIOD) 

FIGURE 6. LUNAR SURFACE TEMPERATURE DATA OF 
LOW, SINTON, MURRAY, AND WILDEY 
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TIME IN DEGREES 

FIGURE 7. FOURIER SERIES APPROXIMATION TO SURFACE TEMPERATURES 



FIGURE 8 





LOSS TANGENT 



THERMAL DIFFUSIVITY (cm 2 /sec) 


FIGURE 9. LOSS TANGENT VERSUS THERMAL DIFFUSIVITY; 
LABORATORY MEASUREMENTS COMPARED TO CALCULATED 
VALUES BASED ON INFRARED DATA, SET I 
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LOSS TANGENT 



THERMAL DIFFUSIVITY (cm 2 /sec) 


FIGURE 10. LOSS TANGENT VERSUS THERMAL DIFFUSIVITY; 
LABORATORY MEASUREMENTS COMPARED TO CALCULATED 
BASED ON INFRARED DATA, SET II 
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FIGURE 11. DIELECTRIC CONSTANT VALUES BASED ON 
INFRARED DATA, SET I 



FIGURE 12. DIELECTRIC CONSTANT VALUES BASED ON 
INFRARED DATA, SET H 


25 





2 4 

O 

O 

E 3 

i— 

o 

LU 

-j o 
uj 2 

Q 


BASED ON INFRARED DATA SET 33 


is 


BASALT POWDER , X = M» cm 


r* PUMICE POWDER 

\ X = 1.1* e» 

k\ 


? * ! 


co ~ 

g§ 

CO ' 

o 11 

-I -< 


BASED ON INFRARED DATA SET I 


FIGURE 13. DIELECTRIC CONSTANT VALUES, X = o. 80 cm 



FIGURE 14. COMPARISON OF LUNAR MICROWAVE BRIGHTNESS 
TEMPERATURES WITH CALCULATED SUBSURFACE 
TEMPERATURES BASED ON INFRARED DATA, SET I 
(SINTON, LOW, AND SAARI) 
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FIGURE 14. (Continued) 



FIGURE 14. ( Continued) 
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TEMPERATURE (°K) TEMPERATURE (°K) 



FIGURE 14. (Continued) 



FIGURE 14. (Continued) 
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FIGURE 14, (Continued) 






TEMPERATURE (°K) TEMPERATURE (°K) 




FIGURE 14. (Concluded) 
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FIGURE 15. COMPARISON OF LUNAR MICROWAVE BRIGHTNESS 
TEMPERATURES WITH CALCULATED SUBSURFACE 
TEMPERATURES BASED ON INFRARED DATA, SET H 
(SINT ON, LOW, MURRAY, AND WILDE Y) 



FIGURE 15. (Continued) 
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FIGURE 15. (Continued) 
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FIGURE 15. (Continued) 
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